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We consider the problem of late-time isotropization in spatially homogeneous but anisotropic 
cosmological models when the source of the gravitational field consists of two non-interacting perfect 
fluids - one tilted and one non-tilted. In particular, we study irrotational Bianchi type V models. By 
introducing appropriate dimensionless variables, a full global understanding of the state space of the 
- . . gravitational field equations becomes possible. The issue of isotropization can then be addressed in 

^— , a simple fashion. We also discuss implications for the cosmic "no-hair" theorem for Bianchi models 

' when part of the source is a tilted fluid. 

O 
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O ■ I. INTRODUCTION 

<N . 

' Since the Bianchi models were introduced into cosmology they have been the most studied generalizations 

i of the spatially homogeneous and isotropic Friedmann-Lemaitre (FL) models. The Bianchi models are spatially 
homogeneous but anisotropic, and thus well suited for studying the effects of anisotropic expansion on the evolution 
\^ of the universe. The complexity of a Bianchi model is determined by its three-dimensional symmetry group G3, 
in conjunction with the chosen matter model. Since the present universe seems to be very well described by a FL 
model, our main interest is models that can be "close to a FL model" at late times. In this paper we will restrict our 
: consideration to models that may approach the open FL model at late times. Within the class of Bianchi models, 
only models of Bianchi type V and Vllh can possibly have this behavior (see Collins & Hawking 1^). This follows 
since the open FL metric admits a G3 of these particular Bianchi types. Here we will focus on the type V models. 
^ The notion of a model approaching a FL model at late times is often referred to as late-time isotropization. 
^ \ As regards the matter description, most studies of Bianchi cosmologies use non-tilted barotropic perfect fluids with 
. a linear equation of state. The book by Wainwright & Ellis Q is basically devoted to such models. In the so-called 
Pf*' non-tilted models, the fluid 4- velocity is orthogonal to the orbits of the isometry group G3. Models in which the 
J> , 4-velocity of the fluid is not orthogonal to the group orbits are referred to as tilted, and were introduced by King & 
' Ellis Q. These models have been studied by, for example, Collins & Ellis ||], Hewitt & Wainwright and Harnett 
. Since the non-tilted models are a subset of the tilted models, the latter can be viewed as simple generalizations 
of the former. Another generalization of non-tilted models is achieved by allowing the source of the gravitational field 
to be a combination of two non-interacting, non-tilted perfect fluids. Since these models can take into account both a 
radiation-dominated as well as a matter-dominated epoch of the universe, they may be considered as more physically 
relevant than single-fluid models. The qualitative behavior of two-fluid models were studied by Coley & Wainwright 
1^ who showed that the models are dominated at early times by the stiffer fluid, and at late times by the softer fluid. 
This is the expected behavior of a universe filled with radiation and dust. 

The next step is to allow one of the fluids in a two-fluid model to be tilted, and subsequently, to let both fluids 
be tilted. Since the dynamics of models containing tilted fluids is more difficult to analyze than that of non-tilted 
models, we will consider the combination of one non-tilted and one tilted fluid. The non-tilted fluid is by definition 
irrotational, but tilted fluids can rotate. We chose, however, to consider the subset of irrotational tilted fluids (For 
a discussion of general single-fluid Bianchi type V models, see Harnett 0). We are thus focusing on the subclass of 
irrotational Bianchi type V models. Both fluids are assumed to satisfy linear barotropic equations of state. 
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Po = (7o - 1)Mo , (1.1a) 
Pt = (7t - l)Mt , (1.1b) 

where po,t are the pressures of the fluids, Ho,t the energy densities, and 7o,t are constants that satisfy < 7o,t < 2 
with 7o ^ 7t. We exclude, however, the specific values 2/3 and 2 since these models behave qualitatively different. 
The indices o and t refer to the orthogonal and the tilted fluid respectively and will be used throughout the paper. 
We will also assume that the energy densities are non-negative, i.e., 

Mo > , (1.2a) 
> . (1.2b) 

To discuss isotropization in detail, we need a well-defined notion of when a model is "close to a FL model" . For 
models with a single non-tilted fluid, a vanishing fluid shear defines the FL models (see, for example, section 2.4 in [Q). 
It is, however, not sufficient to demand that the fluid shear itself should approach zero at late times, since this will 
occur in any single- fluid Bianchi model, irrespectively of whether the model isotropizes or not. As realized already by 
Kristian & Sachs |^ , the appropriate quantity to consider is the dimensionless ratio formed by normalizing the fluid 
shear with the fluid expansion, the so-called dimensionless shear. It measures the dynamical importance of the shear 
compared to the expansion of the fluid. Since there are now two fluids present, the notion of isotropization at late 
times needs to be generalized as follows. We say that a two-fluid model isotropizes at late times if the dimensionless 
shear of both fluids vanish in this limit. This issue was partially addressed by Goliath & Ellis [|lO|, who considered 
models with a tilted fluid and a non-zero cosmological constant. Such models are contained within the models studied 
in this paper if we set 7o equal to zero. We will comment on aspects of the behavior of these models that were not 
addressed in JTo| . In general, vanishing dimensionless shear is not sufficient for a model to isotropize, since models can 
be Weyl-dominated in the future, as is the case for Bianchi type VIIq and Bianchi type VIII models However, 
for Bianchi type V this is not the case . 

The plan of the paper is as follows. In Sec. 0, the field equations are rewritten as a first-order system of autonomous 
ordinar y di fferential equations. Reduced dimensionless variables are then introduced leading to a compact state space. 
In Sec. lH we perform a local analysis of the system of equations. In Sec. Bianchi type I and type V models with 
two orthogonal fluids are studied, while the locally rotationally symmetric type V models are considered in Sec. 
The question of isotropization is discussed in detail. We end with a discussion in Sec. 0. In App. @, the kinematical 



properties of both fluids are given. The relationship to the variables used in [g| is discussed in App. We use units 
such that c = SttG = 1. Orthonormal frame indices are denoted by latin letters a,b, ... . 

II. THE GRAVITATIONAL FIELD EQUATIONS 

The line element for the irrotational Bianchi type V models can be written 

ds^ = -~dt^ + Di{tfdx^ + e-^"" [D2{tfdy^ + Dj,{tfdz^] , (2.1) 

(see, for example, Ellis & MacCallum jl^). In the models we consider there are two preferred timelike congruences. 
First we have the congruence associated with the normal to the spatial symmetry su rfac es, which is also, by definition. 



the congruence of the 4- velocity of the non-tilted fluid, u". The time variable t in (2.1) is chosen so that the normal 
to the symmetry surf aces is The other preferred congruence is that of the tilted fluid. The specific form of the 



line element, Eq. (2.1), guarantees that this fluid is irrotational. Hence, the 4-velocity of this fluid, m", is constrained 

by the field equations to be of a particular form, and can conveniently be parameterized in terms of the so-called tilt 
variable v according to 

< = r(l,t;,0,0) , r = (1 - i;2)-i/2 , (2.2) 

where the value v = corresponds to a non-tilted fluid. The energy-momentum tensor of each fluid is 

= + (7o - 1)?7"1 Mo , (2.3a) 

T,^" = [j.utui + (7t - 1)7^^"] Mt . (2.3b) 

The assumption that the two fluids are non-interacting leads to equations of motion of the form 

VaT^'' = = VaT^^' ■ (2.4) 
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We choose to parameterize the gravitational field using the expansion, 6*, and the non-vanishing components of the 
shear-tensor a±, of the normal congruence (and thus of the non-tilted fluid). They are given by 



-\n{DiD2D3) , C7+ 



2dt 



D2D3 



2 dt 



(2.5) 



The expansion and shear of the tilted fluid can be written as functions of these variables, see App. A. The gravitational 
field equations. 



Gab — T 



ab 



Tl 



ab 



and the equations of motion of the fluids, Eq. (2.4), become 

Evolution equations 



V 



l/)2 2/2,^2' 



-{9-2vBi)a+ , 
-i(0-2a+)Bi , 

V{1 - w2) 



- 2(370 - 2)Mo - 2 1 + (^, _ 1),2 



3[l-(7t-l)^;2] 



[2a+ + (37t - 4)0 - 6(7t - l)Biv] , 



(2.6) 



(2.7a) 

(2.7b) 
(2.7c) 
(2.7d) 

(2.7e) 
(2.7f) 







Constraint equation 

7tw/i„ + 2(1 + (7t - l)v^)Bia+ . 



(2.7g) 



where we have introduced 



Defining equation for fin 

3fin - (7+2 - cr_2 _ g^2 _ 3^^ ^ 



Bl - i?r' : 

1 + (7t - 

= n , 



(2.7h) 

(2.8a) 
(2.8b) 



and a dot denotes differe ntiat ion with respect to t. 

From the form of Ec^. (2.7h), it is now clear why our choice of variables is a good one. The assumption in Eq. (|l.2|), 
in conjunction with Eq. ( |2.7h| ) and Eq. (2.8b), implies that 6 cannot change sign. Theref ore w e can, without loss of 
generality, assume that > 0, i.e., we restrict ourselves to expanding models. From Eq. ( ^.7h ) it also follows that 9 
is a domin ant quantity. We then introduce bounded dimensionless "^-normalized" variables for which the system of 
equations (2.7), is reduced as far as possible, and for which the state space is compact. These variables are defined by 



^^--f' ^""T' ^^"-"0^- 

The subsequent introduction of a dimensionless time variable r, which satisfies 

dr _9 
'dt ~3 ' 

leads to a decoupling of the 0-equation, which can be written on the form 

3(2 - 7t) + (57t - 6)«2 



2_2A2_3(2-7o) 



2 [1 + (7t - l)v^ 



(2.9) 



(2.10) 



(2.11) 



where the prime denotes differentiation with respect to t. The parameter q is the deceleration parameter associated 
with the normal congruence and the non-tilted fluid. The remaining equations can now be written on dimensionless 
form as follows. 
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Evolution equations 



A' 



n' 



-(2-9)S_ , 
(g + 2S+)A 



1- (7t-lK 
[2q - (37o - 2)] 



[2E+ + 37t-4-2(7t-l)Ai;] , 



(2.12a) 
(2.12b) 
(2.12c) 

(2.12d) 

(2.12e) 



Constraint equation 

= jtvn,, + 2[i + (7t - i)v^]aj:+ 



(2.121) 



Defining equation for 

1 - Ei - - ^2 _ 



(2.12g) 



The set of equations ( 2.12| ) shows that the irrotational Bianchi type V models with one orthogonal and one tilted 
fluid is governed by a system of five autonomou s ord inary differential equations subject to one constraint. The 
dimension of the state space is thus four. The set (2.12) is invariant under the discrete transformations 



(2.13) 



so we can, without loss of generality, restrict ourselves to the invariant set defined by A > and S_ > 0. The 
varia bles £ +, S--, A, v, and flo therefore satisfy < E^, A, w^, fto < 1. There is a number of important subsets of 
Eqs. (|2.12D: 



(i) The three-dimensional subset defined by u = 0, which describes a Bianchi type V universe with two non-tilted 
fluids. This case is thus included in the general study of Coley & Wainwright We will, however, consider 
these models in Sec. IV for the purpose of comparison with the models where one fluid is tilted. The subset 

V = also contains Bianchi type I models and the open FL model as special cases. 

(ii) The three-dimensional subset — 0, which corresponds to locally rotationally symmetric (LRS) models (see, 
for example, |13|). This subset turns out to be very important for the evolution of the general irrotational type 

V models and will be considered in detail in Sec. ^ 

The next step in the analysis is to consider the equilibrium points of Eq. ( ^.12| ) . This is done in the next section. 



III. QUALITATIVE ANALYSIS 



In Table |, we present the equilibrium points of the system (2.12). The corresponding eigenvalues are given in Table 
||. Since the constraint, Eq. ( 2.121 ), cannot be solved globally in an analytic way, it will be treated locally (see, for 
example, Q). 

Some of the equilibrium points correspond to exact solutions of the field equations. For example, the equilibrium 
points and correspond to flat FL models in which the tilted fluid and the non-tilted fluid is dominant, respec- 
tively (The "tilted" fluid is in fact non- tihed at F° since u = 0, but we refer to it as the tilted fluid for simplicity). 
The point M° is the Milne model, while the equilibrium set K'^ corresponds to Kasner-like models. The point M 
corresponds to flat space and coincides with Af" for 7t = 4/3. We also note the appearance of the equilibrium set F^ 
for the specific value 7t = |- It is associated with a line bifurcation that transfers stability between the equilibrium 
points F^ and F^ . The points correspond to particular Kasner solutions. There is also a number of equilibrium 
points for which the tilt is extreme (w^ = !)• Whether these equilibrium points correspond to exact Bianchi soluti ons 
or not seems to be an open question (see comment on p. 4245 in |^). We note that the constraint, Eq. ( 2.12i ), is 
degenerate (VG = 0) at the point allowing all five eigenvector directions to be physical at this point. 
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Restrictions 
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1 




M" 








1 










K° 




E^. + e2 = 1 













M 








1 


_ 37t-4 
^ 2(7t-l) 





f < 7t < 2 


C± 


E+ = -i(37t-4) 


E_ = ±iV3(2-7t)(37t-2) 





V 







F^ 











V 


1 


7t = 1 


-'■ o 








±1 


1 






M± 








1 


±1 











E2. + e2 = 1 





±1 







n 


T.+ = -l + A 





A 


1 








TABLE I. Equilibrium points for the irrotational Bianchi type V models with one orthogonal and one tilted perfect fluid. 
Unless indicated, all equilibrium points and sets are in the physical part of state space for < 7o,t < 2, 7o,t ^ 2/3. 





Eigenvalues 


EUm. 


F^' 


-1(2-70 
-1(2 -7o) 


-1(2-70 
-1(2 -7o) 


4(37t-2) 
i(37o-2) 


3(7t - 7o) 




V 


fS 


37t -4 


-3(7t - 7o) 




M° 


-2 


-(37t-2) 


37t -4 


-(37o - 2) 






K° 





3(2 - 7t) 


2E+ + 37t - 4 


3(2 - 7o) 




A 


M 


-2 


2-7t 
7t-l 


(37t-4)(57t-6) 
(7t-l)(97t-10) 


-(37o - 2) 






C± 


3(2 - 7t) 








3(2 - 7o) 




A 


pv 
^ o 


-1(2 -7o) 


-1(2 -7o) 


5(370-2) 









' o 


-1(2 -7o) 


-1(2 -7o) 


|(37o-2) 


2{37t-4) 
2-7t 






M+ 


-2 





2 


-(37o - 2) 




s+ 


M- 


-2 


-4 


2(57t-6) 
2-7t 


-(37o - 2) 




s+ 




2(1 + E+) 





-2^[2S+ + (37t-4)] 


3(2 - 7o) 




A 


H 





-2(1 + E+) 


2(1 - 2S+) 


-(4E+ +370 -2) 




A 



TABLE IL Eigenvalues for the equilibrium points of the irrotational Bianchi V models, with one orthogonal and one tilted 
perfect fluid. 
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IV. TWO ORTHOGONAL FLUIDS 



Setting u = in the constraint Eq. ( 2.12j ), implies either A = 0, S+ 7^ (Bianchi type I models) or — 0, A ^ 
(Bianchi type V models). Without loss of generality we can assume that -ft > Jo- The boundary subsets of the state 
space for these two classes of models are given by the two invariant submanifolds f2o = and lln = 0, which describe 
the corresponding one-fluid models. The dynamics of the Bianchi I state space is shown in Fig. |l|, while the dynamics 
of the type V models is shown in Fig. 

For type I models there are two sources. The equilibrium point F° gives rise to a single orbit ending at F^. It 
corresponds to a flat Friedmann model with two orthogonal fluids. The other source is the equilibrium set of 
which each point is associated with a one-parameter set of orbits. Therefore this equilibrium set describes th e ge neric 



behavior at early times. The future attractor of these orbits is the point Thus, from Eqs. ( [A^ ) and (A9), all 
orthogonal two-fluid Bianchi type I models isotropize. 

For type V models there are three sources. The equilibrium point is associated with a one-parameter set of 
orbits (characterized by S_ = 0), which corresponds to open Friedmann models with two orthogonal fluids. They all 
end at the equilibrium point M°. The other two sources are the two equilibrium points belonging to the set K'^ with 
E_ = ±1. Both of these are associated with two-parameter sets of orbits, and thus describe the generic early-time 



behavior. They are future attracted to A/°. Consequently, from Eqs. ( [A^ ) and (A9), all orthogonal two-fluid Bianchi 
type V models isotropize. 

From the analysis of the orthogonal two-fluid models, it is clear that the generic behavior is described by equilibrium 
points associated with two-parameter sets. Thus, in what follows, we will focus on such equilibrium points. Note that 
the variable x oi Coley & Wainwright |^ (see App. ^ is a monotone function when both the fluids are orthogonal. 
This is no longer the case when one of the fluids is tilted. 



V. BIANCHI TYPE V LRS MODELS WITH ONE TILTED AND ONE ORTHOGONAL FLUID 



As Bianchi type I mod els do not allow the combination of one non-tilted and one tilted fluid as a source due t o 



the constraint Eq. ( 2.12] ), we focus on the Bianchi type V models. Since J7o,rin > imply q < 2 by Eq. ( 2.11 ), 
the evolution equation for E_ implies that S- is a monotone decreasing function along all orbits with I]_ > and 
riojf^n > 0. This fact significantly restricts the evolution at late times. It implies that limT-_^oo S_ = 0, for all 
orbits with r2o,^^n > 0. This can be proven along the lines used to prove the similar statement for tilted single-fluid 
models, see |^]. The asymptotic behavior as r — > 00 is thus contained in the three-dimensional invariant set E_ — 
corresponding to the LRS models. 

The three Kasner circles K'^ and each reduce to two equilibrium points in the LRS submanifold, namely points 
for which S+ = ±1. We denote these equilibrium points K'^, f^nd K^, where the subscript distinguishes between 
the two signs of S+. The equilibrium points for which v and S+ have the same sign (collectively denoted K^) are 
not located on the boundary of the interior of the LRS submanifold. Consequently, we do not need to consider them 
when studying the dynamics. 



The state space for the LRS submanifold with 2/3 < 70 < 7t < 2 is presented in Figs. |3|-|^. The effect of changing 
the equation-of-state parameters so that 70 > jt is that the flow along the orbit between F^ and F° is reversed. 
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FIG. 3. The state space of Bianchi type V LRS models 
when with 2/3 < 70 < 7t < 6/5. 




FIG. 5. The state space of Bianchi type V LRS models 
when 7t = 4/3 and 2/3 < 70 < 4/3 



FIG. 4. The state space of Bianchi type V LRS models 
submanifold when 7t > 70 with 2/3 < 70 < 4/3 and 
6/5 < 7t < 4/3. 




K f" k; 

FIG. 6. The state space of Bianchi type V LRS models 
when 7t > 70 with 2/3 < 70 < 2 and 4/3 < 7t < 2. 



Similarly, < 70 < 2/3 results in a stability change along the orbits from the equilibrimn points F* to M*, where 
* G {0, +, — }. Note that the special case 70 = corresponds to a cosmological constant. Consequently, the LRS state 
space for < 70 < 2/3 is contained in Figs. 9-11 of |jl^. Finally, if < 7t < 2/3 the flow changes along the orbits 
F°-M° and Kl~K:^. 

The sources and sinks for var ious equat ions of state are summarized in Table The isotropization properties, 
which can be found from Eqs. ( [A^ ) and ( |A6| ), are also listed. For 7t < 6/5, all models isotropize. for 70 > 2/3, 
6/5 < 7t < 4/3 there is a class of solutions of non-zero measure that does not isotropize, namely the class of solutions 
associated with M~. For 70 > 2/3, 7t = 4/3, all models isotropize. For 7t > 4/3 no models isotropize. Note that 
the physically interesting combination of one dust fluid and one radiation fluid always isotropizes, regardless of which 
fluid is tilted. 

When 7o < 2/3, 7t > 4/3, the future attractors are F^. Solutions corresponding to orbits approaching these 
equilibrium points do not isotropize, see Table II], even though the final states are inflationary in the sense that the 
deceleration parameter associated with the normal congruence q = 5(370 — 2) is negative. Thus, in some sense it 
seems to be misleading to refer to solutions corresponding to orbits approaching these equilibrium points as "asymp- 
totically Friedmann models" , although the solutions corresponding to the points themselves are Friedmann models. 
In particular, for the case of a cosmological constant 70 = 0, the solutions corresponding to F^ are de Sitter models 
pO| . This is consistent with the cosmic "no-hair" theorem for Bianchi models [H. However, the theorem does not 
guarantee that the tilt tends to zero, as pointed out by Raychaudhuri & Modak~|[l5|| . From our analysis, it is clear 
that this cautionary note is crucial for these models. As the tilted fluid does not become orthogonal at late times, 
the expansion-normalized shear of the tilted fluid does not vanish at late times. This stresses the point made in |l^] 
that one should be cautious about the isotropization of tilted models. As seen in Table [II, this is in fact the generic 
behavior for models with < 70 < 2/3 and 7t > 4/3. 

For general irrotational Bianchi type V models, there is no a priori reason that St_ — > w hen I]t+ ^ 0. However, 
this is indeed the case for the equilibrium points in question, as can be seen from Eq. (A8), noting that S_ 0. 
Thus, the analysis of isotropization for the LRS submanifold holds for the general class of models as well. 



VI. DISCUSSION 



In this paper we have continued the study of irrotational Bianchi type V cosmologies, using the dynamical systems 
approach initiated by Hewitt & Wainwright ||] and Coley & Wainwright The source of the gravitational fleld 
has been taken to be two non-interacting fluids, one orthogonal and one tilted. Such models can describe a universe 
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sources 


< 7t < 2/3 
2/3 < 7t < 2 


kI, n 




Sinks 

< 7o < 2/3 2/3 < 7o < 2 


< 7t < 2/3 

2/3 < 7t < 6/5 
6/5 < 7t < 4/3 
7t = 4/3 
4/3 < 7t < 2 


-ft" (7t < 7o) all isotropizc 
(7t > 7o) 

F° all isotropize 
Fo all isotropize 

none isotropize (except for t; = 0) 
none isotropize 


-F" all isotropize 

all isotropize 
M°, M" M° isotropize 

M°, M" all isotropize 
M, M~ none isotropize 



TABLE III. Sources and sinks for the Bianchi type V LRS models with two fluids. 



where one of the fluids models the contribution of radiation to the energy density of the universe, and the other the 
matter content. 

We have found that, although the orthogonal two-fluid models isotropize, this is not necessarily the case when one 
of the fluids is tilted. Thus, depending on the equation-of-state parameters 70 and "ft, it is possible to find cases 
for which all or a subset of the solutions are anisotropic to the future. In particular, there are models which are 
inflationary, but do not isotropize. However, we emphasize that the cases of dust plus radiation always isotropize. 

It should be stressed that the relevant quantities when determining whether the shear dies away are shear quantities 
normalized with respect to the expansion associated with each fluid. 
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APPENDIX A: FLUID PROPERTIES 



Here we present the kinematical properties of the two fluids in terms of the variables used to parameterize the 
gravitational fleld. For the non-tilted fluid we have 



Fluid expansion 



Fluid shecir 



(To = <J+ + (T- 



(Al) 
(A2) 



The dimensionless shear E„ of the non-tilted fluid is thus 



V2 _ ° _ v2 I v2 



The fluid properties of the tilted fluid are as follows. 
Fluid expansion 

'W-QvBi+v^{2a+-9] 



1 



3[l-(7t-l)t^2] 



Fluid shear 



l-t;2 



2vv 



3- 2wA + t>2(2i;+ - 1) 
3[l-(7t-l)t'2] 



+ - Biv{2a+ - Blv) - ^(c7+ - B-^v) + 

1 — 



{l-v 



2\2 



(A3) 



(A4) 



(A5) 
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The dimensionless shear St of the tilted fluid can be written 



= ^ - + S?_ , (A6) 



01 



where 



y _ 1 I 3(l + E+-.A)[l-(7t-iy] 
_ 3[l-(7t^l)^^^]S- 

" 3 - 2^;^ + (2S+ - ' ^^^^ 



(see Eqs. (A3), (A4) and (A6) in g). Note that when both fluids are orthogonal (u = 0), Eq. (A6) reduces to 

Y,l = Yi\ + Y?_. (A9) 

APPENDIX B: THE VARIABLE x 

In their study of two orthogonal fluids, Coley & Wainwright introduced the following variable instead of VLq: 

The evolution equation for % becomes 

X' = - ^l\^^^^l)^2 {37o [1 + (7t - 1)^'] - 7t [3 + - 2t;(^ + z;S+)] } . (B2) 
For the submanifold v — Q corresponding to two orthogonal fluids, the above equation simplifies to 

x' = -^(l-x')(7o-7t) . (B3) 
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